General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




MSC 


INTERNAL TECHNICAL 


NOTE 


GENERATION OF RANDOM VECTORS WITH 
KNOWN MEAN AND COVARIANCE 



* 


NATIONAL 



AERONAUTICS AND SPACE ADMINISTRATION 
MANNED SPACECRAFT CENTER 
HOUSTON, TEXAS 



June 6, 


( ) 



MSC INTERNAL NOTE NO. 65-ED-21 


GENERATION OP RANDOM VECTORS WITH 
KNOWN MEAN AND COVARIANCE 


Prepared by: 
i 



'liAL (lAAi 

Fred M. Speed / 
AST, Data Systems 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
MANNED SPACECRAFT CENTER 
HOUSTON, TEXAS 

June 8, 1965 


GENERATION OP RANDOM VECTORS WITH KNOWN MEAN AND COVARIANCE 


By 

' Fred M. Speed 

MANNED SPACECRAFT CENTER 

* 

• SUMMARY • 

In simulation studies, it is necessary to introduce 
normally distributed random "errors” to the data in order 
to simulate the actual conditions. These random "errors” 
are normally distributed random vectors with known mean 
and covariance. The algorithms currently being used (ref. 3) 

t 

to generate these vectors encounter one of the two following 
problems. 

(1) The algorithm can not handle singular covariance 
matrices (a singular covariance matrix occurs 
when the data is completely correlated) . 

(2) The algorithm requires the inversion of a matrix. 

The purpose of this paper is to present an algorithm that 
does not require the inversion of a matrix and can handle both 
singular and non-singular covariance matrices. The proposed 
algorithm should be more efficient than the current algorithms 
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INTRODUCTION 

% 

This paper is divided into four parts. The first part 
gives some definitions and theorems that will be needed in 
the development of the algorithm. The second part contains 
the development of the algorithm. The third part describes 
a computer program that' generates ' the desired random vectors,, 
The last part contains Appendices A and B. Appendix A contains 
a listing of the computer program GERN and presents several 
examples illustrating the use of the program. Appendix B 
contains listings of the subroutines that make up the N(0, 1) 
random number generator packet. 


i 

SYMBOLS 

X 

column vector 

Y 

column vector 

R 

positive semi-definite symmetric 


matrix 

Other capital letters 

matrices, unless otherwise stated 

a t 

transpose of A 

A' 1 

inverse of A 

Y is distributed 

Y is from a normal population with 

according to N ( y , R) 

mean y and covariance matrix R 
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E(Y) 

COV(Y) 

X 

' <P 
D 

DEFINITIONS AND THEOREMS 

Definition 1. Two n x n real matrices A and B are said 
to be congruent if there exists a non-singular matrix P 
such that 

j 

I T 

! PAP = B . 

Since congruence is a special case of equivalence * the 
matrix P can be obtained from elementary row-column 
operations (ref. 1). 

Theorem I. (ref. 1) Every n x n real symmetric matrix 
A of rank r is congruent to a diagonal matrix, whose 
diagonal elements consist of r ones and n-r zeros. 

Theorem'll, (ref. 2) Let X be distributed according to 
N(<£, H) . If Y = A.X, then Y is distributed according to 
N( <f>, AHA T ) . 


expected value of Y 
covariance of Y 
identity matrix 
null vector 
diagonal matrix 


Theorem III* (ref. 2) If Y = AX + y, where A is a 
constant matrix and yi is a constant column vector, then 

E(Y) s AE(X) + p 
COV(Y) = GOV (AX) 

The following theorem is an immediate consequence of 
matrix theory. 

Theorem IV. If D is a diagonal matrix, whose element" 
consists of zeros and ones, then D = D’D. 

! 

• THE DEVELOPMENT OF THE ALGORITHM 

Suppose it is necessary to generate an n x 1 random vector 
Y with mean <j> and covariance R. The following theorem 
provides a new method to obtain this Y. 

Theorem V. Let X, an n x 1 vector, be distributed according 
to N ( 4» , I). If R is a n x n given covariance matrix, then 
there exists an n x n matrix A, such that if Y = AX, then Y 
is distributed according to N(<{,, R) 

Proof: Let the rank of R be r < n. By definition, R is 

a real symmetric matrix. Thus, by Theorem 1, there exists 

T 

a non-singular matrix P such that PRP ~ D, where D is a 


diagonal matrix whose diagonal elements consists of r ones 
and n-r zeros. 


Theorem XI implies that a necessary and sufficient 

condition that Y - AX be N( , R) is that AA T » R. It will 

T T 

be shown that A = RP is one matrix such that AA = R. 

Since PRP T = D, then 

(i) RP T = P -1 D 

rn 

(ii) PR = D (P 1 ) 

(ill) R = P" 1 D (P T ) 

Thus, if A = RP T , then 


m T T 
AA 1 = RP PR 1 


T 

s RP PR 


= P" 1 D D(P T ) 


-1 


= P*" 1 d(P T ) 


= R 

Thus, if A - RP T , then Y - AX = RP T X is R) and the 

theorem is proved. 

It has been shown (ref. 1) that P can be found by forming 
the augmented matrix [R, I] and operating on this matrix 


I 



with elementary row-column operations in such a way that 

R is diagonalized. The result of these operations on 

* 

[R, X] is [D, P], Since X is distributed according to 
W(<j> a I)^ the elements of X are readily obtainable from a 
N(0, 1) random number generator. Thus,. given the covariance 
matrix R, the random vector Y can be obtained rather easily. 


SUBROUTINE GERN 



GERN is a FORTRAN IV subroutine used to generate random 
vectors having a known mean and known covariance matrix. 
All computations are done in singlet-precision floating 
point arithmetic. The method used is described in the 
previous section. 


Calling Sequence 


Call GERN (N, R, Y, A) 
where : 


N 


size of R 


O 


R 


covariance matrix. R is 
dimensioned R(50, 50) 
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- ♦ «ir% 4 »in » trrrt e i’L. XKjsi,, -i,,^. 


Y n x 1 random vector from a N($, R) 

*> 

population, Y is dimensioned 
YC50) 

A a matrix such that Y - AX is 

distributed according to N(f, R) , 

A is dimensioned A(50, 50). 


Error Messages 


R is not positive semi-definite, a negative number will 
occur on the diagonal of R during the row-column operations. 
When this occurs, the following message is printed: 

"THE ORIGINAL MATRIX IS NOT POSITIVE SEMI -DEFINITE . 11 


Restrictions 


N < 50 


N(0, 1) Random Number Generator 


See Appendix B 


Method 


Given : 


R 


y 


Construct: [R, I] 
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Operate on [R, I] with 'elementary rovr-column operations 
to obtain *[D, P] . 

Construct : X ™' 'M ^ a j 33? 

Compute : A = RP T 

Compute; (Y)~ AX 

y-v u(o 

CONCLUSION 

Theorem V provides a very simple method to obtain random 
vectors that can be used to simulate observations that are 
highly or even completely correlated as well as observations 
that are independent . 



REFERENCES 


1. Ayreb Jr., Frank: Theory and Problems of Matrices. 

*• 

New York, Schcrem Publishing Co., pp. 40~*12, p. 115. 

2. Anderson, T. VI.: An Introduction to Multivariate 

Statistical Analysis, New York, John Waley and Sons, 

Inc., p, 2,5. 

3 . Tapley, b, D. : Odell, P. L. : A Study of Optimum 

Methods for Determining and Predicting Space Vehicle 
Trajectories and Control Programs. Contract NAS 9-2619. 

b. Brunk, K. D. : An Introduction to Mathematical Statistics. 

Boston, Gin and Co., pp . 86-90. 


APPENDIX A 
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This appendix contains three examples and a listing of 
the Subroutine GERN. 

Example 1: Suppose it is necessary to simulate a 5 x 1 

radar observation vector Y } where 


Y 


range 

azimuth 

time 

elevation 

range-rate 


In order to simulate the actual conditions , random 
"observation errors" must be added to the vector Y. Let 
these "errors" be from a N(<f>, R) population 3 where 


1.0000 

0.5576 

1 ' 

0 .4641 

0.8197 

0.2333 

0.5576 

2.0000 

0.1719 

0.2516 

0.2265 

0.46^1 

0.1719 

3.0000 

0.0264 

0.0334 • 

0.8X97 

o . 2516 

0.0264 

4.0000 

0.9608 

.0.2333 

0.2265 

0.0334 

0.9608 

5.0000 - 


Subroutine GERN computes the matrix 


n 



'i.oooo 

f 

0.0000 

0.0000 

0.0000 

0.0000 ' 


0.5576 

1.2996 

0.0000 

0.0000 

0.0000 

A = 

0 . J| 6 i) X 

-0.0668 

1.6673 

0.0000 

0.0000 


0.8197 

-0.1581 

-0.2186 

1.8042 

0.0000 


. 0.2833 

0.0742 

- 0,0419 

0.4279 

2.1806 _ 


This matrix is then multiplied by X, where X is 
distributed according to N(<f>, I). The result of this 
multiplication is the random "error 11 vector 



”- 1 , 0675 ’ 
1.0053 
-I.7816 
. - 1.6026 
_~ 4 . 24 l 7 _ 


The vector Z is added to the vector Y to obtain -the' 
simulated radar observation vector, 


Example 2 : Suppose> in Example 1 , it is assumed, th 

can be measured without error , ....This 
the vari an 

Under this assumption, the . 






1 


from a N( <f> > R) population, where R might be 


R = 


1.0000 

0.2248 

0.0000 

0.9471 

0.4625 

0.221)8 

2.0000 

0.0000 

0.0865 

0.6449 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.9^71 

0.0865 

0.0000 

4.0000 

0.2663 

i 

0.4625 

0.6449 

0.0000 

0.2663 

5.0000, 

iatrix A 

is computed 

to be 



1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.2248 

1.3962 

0..0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.9^71 

-0.0905 

0.0000 

1.7591 

0,0000 

0.4625 

0.3873 

0.0000 

- 0.0777 

2.1517 


and the random "error” vector Z is 


0.5879 

0.2784 

0.0000 

1.7551 

•0.5159 
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Example 3: Suppose it is necessary to simulate a 6 x 1 

observation vector Y, where 


Y « 


*1 

y 2 


y 4 

*5 

H 


and where yg = l y 


1=1 


This is a case where an element of the observation 
vector is completely correlated with the other elements. 


the "errors” 

be from 

a NU, R) 

population where 



2.000 

0 . if 11 

1.334 

- 0.097 

1.612 

5.259 


0 All 

O 

o 

o 

• 

-=x 

- 0.238 

-0.684 

-0.656 

2.832 


1-334 

-0.238 

6.000 

-1.590 

1.024 

6.530 

R = 

-0.097 

— 0 . 6 8 if 

- 1.590 

8.000 

-1.226 

4.401 


1.612 

- 0 . 656 ’ 

1.024 

-1.226 

10.000 

10.755 


5.259 

2.832 

6.530 

4.401 

10.755 

29.779 


» 


Note that the last row of H is the sum of the first 
five rows of R and hence R is a singular covariance matrix 

Subroutine GERN computes the matrix 


A = U 


"l.ltlil 

0.000 

0.000 

0.000 

0.000 

0.002 

0.290 

1.978 

0.000 

0.000 

0.000 

0.002 

0.91)3 

-0.258 

2.2115 

0.000 

0.000 

-=r 

O 

O 

O 

-0.069 

-0.335 

-0.718 

2.711) 

0.000 

0.001 

1 • Xlit 

-0.1198 

-0.080 

-0.506 

2.861 

0.000 

3.719 

0.886 

1. 1)H7 

2.208 

2.861 

0.005 


The random ’’error” vector Z is calculated as before and 
the result is 

- 1.581 
-0.051 
3.311 
Z = -0.496 

.688 

1.874 

A listing of Subroutine GERN follows. 
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tlttFTw G r WN 

S l • n 0 1 • T T M r Gfr °N( N*k» Y # M) 

c ^ AS MY p • '*• SPEED MSf* • ApRjt #19*S 

C CALL IMG SFVJfcMCE 

C CALL GF.3N( Ntr<» Y* M 

c fc^ERE 

C N IS IHF SI«' r 9 C .( 

C ft ts HE* COu Ahl AMC“ maTrTx, 9 IS n i mfvST VF.D PIS'SSO) 

C Y IS RAMU^M W£CT0k PPjV A N {/% R ) POPULATION , V IS nTye-NSjON' 

C YCbi!) . 

C A |S A VtAl-tlV S'.ICH IMaT VzAY l c; f^IST 0 ! ^llTED *C r O»DTN^ TO N(/* 

C A IS ’M ME Mbl i) M EO MSOtbfl)* 

DT Mf VS I ON t ( ) * * < bn | 9 h (59 t 1 09 ) » *CR 0 # * 0 ) ( Sn » ?>n ) 

A t'MsrwbO) »H(s°»‘>n> 

CALL 9ESIMIT) 

L = ?*N 
LI = MO 

C CONSTRUCT C«*i) 

00 1 I = t t M 

DO M J=Ll#L 
K = J-M 

A ( 1 , J) = n.o 
A ( K , J ) r 1 • 0 

1 comtimje 

DO 33 1 = 1 • .m 

no 33 J = i » n 
A ( I , J ) z 9 ( I , I ) 

33 COMflMJE 
Ml z M-1 

C COMPUTE (0»P) 

00 * T z 1 rMi 
I c (A(T#l))4#t>*U 
0 CONTI MtJ" 

Lu z T ♦ 1 

no io < =l 


jn 

CON I 1 MU" 
CO TO 4 


1 1 

COMT 1 VU: 

DO i? Jz 1 » u 



A ( I , J ) Z A ( I , n 

f A ( < » J ) 

\2 

COMT 1 MJE 

no m j = i » m 



A ( J, n = M J. T ) 

i A ( J t S ) 

13 

f ON T I MUE 



4 CONTTMij- 

L ? - Til 

DO 6 <=L2»\ 

D1 z A ( < r T ) 

00 h J= 1 . L 

A ( < • J ) = M<» •) - ( 01 *a( 1 *J))/A(I#I) 
6 COMT IMU- 

no j = l?» m 
a ( i , j ) = n.o 
<!\ CONTINUE 
2 CONTI MJE 
no 7 b i - l t \ 
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reproducibility of the original page is POO 


o 

ic(<v n,tn 

7'J W» 11 f|f IH 14 ) 

J‘»m Cy^v/l(4rtH |Mfc O^oTmAl vAi^lX T S NOT PO<;TiTv r c;pm \ -n^F j N* U ) 
GO ID IS 
7 H CONl(MJ" 

t)0 ot 1 = ] , M 

DO n 1 J - i # N 




C 

c 



c 


H ( 1 # J > = n.o 
DA hM S in 
N .At> S S ♦ M 

(?2 = 'MI# I) ♦ K< l # K ) *M J*N35) 

bl CONI I MU? 

DO '/H I = 1 1 \ 

no l r SO" MriU # T ) ) 

DO J = J # l 

d-1 A ( I , J ) r A(i f j)/0?i 

coNriM* 6 : 

PI = HAMS^O^l J r ° 

COVPUTE • K°1 

no *1 Ir 1 #.\! 

f)0 f \\ J=l»\ 

H ( I » J ) - 0 . U 
DO \? \:1 

N3S r N ♦ K 

*2 H(l,J) = MU, I) ♦ K( | r * ) ♦ A ( J,M3S) 
»<t CA'invil- 
COmPI.H F Y =l*PTMY 

DO r,U I = 1 # J 
b'* X C l ) r 3AP,^\j(T) 

DO hS 1 = 1 » N 

Y ( I ) = 0.0 
DO nO J = J # M 

bn Y ( I ) r Y ( I) *• u ( 1 » J) *Yt J) 

o'-, CO K 'UN‘Jf 

IS CONTI MU? 

UPTURN 

FND 


































reproducibility of the originaTTagm^c? 


WCI Tt It . t • t 4 


ctn*t to ci 
roc 

l>l 


C (I .iMB <1 .J)«c (I ,M •* (i.NJJl 


• C*l A t to 26 

re* 

Ml . 1*1 1 


C O'! TIM.'. 


\ certit to n\ \ rerett to riA . \ ctrcct to «2 

-) ro« ) rod \ — 

/mi.i«i n/ n/ H ...J /mi, 


II |I)SH II • J I ♦ R !| >K) •* I J (II! ) | 


CirtAI TO 64 

re* 

Ml .1*1 1 


I (IMCARMIIT) 


«CK*T TO t} 

re* 

Ml 


CE6CAT TO Cl 


MI|:rill<ll(],JI«l (II 
























APPENDIX B 


’ N(0, 1) Random Number Generator 

The N(0, 1) random number generator packet consist 
of the following subroutines. 


(1) 

BEGIN 

(ii) 

FIX3 

(iii) 

BARNN 

(iv) 

BARNA 

(v) 

LOOKUP 

(vi) 

BLOCK DATA 

( vii ) 

BARN 


All these subroutines are in FORTRAN IV, except BARN, 
which is in MAP. The first six subroutines in the packet were 
programed by K. Oney of Wolf Research and Development 
Corporation, Houston, Texas. The last subroutine, BARN, is 
a SHARE subroutine. 


A listing of each of the subroutines follows.. 
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T1MFIC L T ‘3 T 

flPROhT 1 'IF UKMU ( I ) 

PI VE M^TOM TrtH L «M30) # T A^cf! 1 < l3h) 
COMMON /^\ r r/ T AML r #1^1.11 
C 

CALL lTv.c; ( A ) 

CALL PI XI ( i , **) 

1 = X T 
WFTUHN 
E.no 


M rIF TC FJXM LIST 

miMHO.lT ] Mf F I K n l.l* J» 1 T*‘7) 

c 

TIM = C'L OA T ( Js ° ) 

T I vi | = TIM * 1 ,f-S 

T I V2 = FLOA 1 ( 1 C 1X C T i vi i ) ) *1 ,* ♦ 

T 1 m 3 = 1 TV - | TVi? 

T 1 V4 r T I M3 ♦ 1. fc *-3 

TIVS = FlO^i ( irix(TiMu) )*l. 0 Ff 3 

TI Vt, = 1 I M3 - TI MS 

T I v7 r 1 I M** i , r -5 ♦ 

^FTIN'J 
r \<n 


tIPFTL L T T 

FIlNCllPN PA-<M;j ( T ) 

XT = RA?<M( t . ) * T 

wPTlF(S»3)T#X| 
if (xt , l t • u.mxr = Ai^ivn 
2 FORMAT (?<bx# Fill, 7 ) ) 

1 F(Xt.M . i . n i so ta i 
CALL L^O* U n ( X 1 , xr,) 
flA 3 NN = *0 
PFTIMN 

l xt = xt - i.n 
call LO 0 <U” ( x T » X r * ) 

MA^NN r MG 
PF TUPN 
c^O 


TI^FTC LOO<u D L T ST 

SU^POllT 1 VJF LO«.)XU y (XfP^lN iV) 

C 

^TMENSTOv) TAPL r ( n n) . TA^LFl ( 135 ) 
COMMON /0 N r T/ I ^HL r » 1 A^l El 
C 

F L A G z 0 • 0 
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PLAG1 = O.o 
IF ( x .11. o.t )*u ? 

if t k .Li. o.ojro to s 
x = 1.0 - x 

FI. 1 = 1.0 
GO To 9 

S 11 = I Fi y c V# ion. ) 

1=11-7 

XO = FLOAT (in / too, 

XI = XO f 0,01 
P r (X-XO) / ( v i - y o ) 

4 PS = P*P 
PC = P^*P 
PF = »C* 3 
PS = Pf * 3 

AV? = n, AAX.iVT*'*** -P ♦ (P*p,ll-P*;fS,0-or#R,0-M'F*c # 0-ns) 

AMI Z -O.U1b6fi*l 7 C -1 i ( *>♦ 1 2 . 0-PS* 1 6. 0 + PC ♦°F*u , n-PS ) 

AO = ), ^333333^.-1 t ( 1i...l-P*4,0- oc ;MS.n*P r *S.O4PF*3. 0-PS) 
At = 0.*i333.J3.V»t -t * (p# 1 2,u + Pc;»0,n-Dc*7.n-PFi''.04PS) 

M> z -0.'4l*bSt>*7r-1 ♦ (F»f,.04Pc;-nr^7.0-pC4 n e>) 

A3 = o. G33333.V»t.-P 4 (p%4 .n»Pfitb. 04PS) 

if (Flag .\' r . 0. 0)00 .TO t 

PhTNv r amo. t>» n l r ( T -P ) 4 A v \ ♦ T ApLF (I - 1 ) ♦ * 0*T*P|F(n 

1 t A ] ♦ 1 AOLt ( I ♦ 1 1 4 4£«1*t'.LF(I+?) ♦ A34TARLF ( 1*3) 

IF ( h L A G 1 ,Fj, o,n>Gj to s 
PHTN\/ r -PHiMv- 

V r i . n - v 

n < F T 1 1 r? X 
? / r alo; : ( x ) 

I x = l f i x ( 7 ♦ 1 u . ) / d 

XO = (PLOAH IX) 4?.o - P.n) ♦ 0,1 

vi z yo 4 o , ? 

P = ( 7-xo) / ( Vl-VO) 

I r lAWS(lX) - 7 
FLAG = 1.0 
GO TO u 

3 P M T N v - ^v?4 T,» r ’l. t 1 ( 1 » ?) -4 AVI ♦TAPL r l ( T 4-1 ) ♦ A0*TARl. F\(T) 4- 

J A 1 ♦! AHLr-1 ( T-l) 4 A?4 TAPL r l ( I-?) 4- A3*TApt.F1 ( T-X) 

I F ( F L A G 1 ,Fj, 0.0)Gp TO 7 
PHI Si v = - 3 H T ,\ v/ 
v r 1 ,0 - V 
7 HTTlK>i 

f NO 


tlPF T (. PAR\A LTST 

F l ) V C T I ONI PAr<\|A f T ) 

XT = H A k M ( 1 , ) 4- T 

IF (XT , G T • 1.0) GO TO 1 



HA^NJA r Xl 
>n TUtfN 

I P.A3NA = XT - 1.0 
HFTU’<\' 

ENO 




iigftODUCmn irv 0 

LOF 7HE original 



U M»»U»C LThT 

‘•l ;'f < 

COVVjom T 

.'‘J MKV|Cf {)\. T*. » t r 
f:AT A ( TA-»l 1 J ( i > 
x -n. i uuu M'9«.r 
x -o , | i*xx \ 7 h*tr 
X -0.O79U l 0*ii. r 
X - 0.91 711 VII 7 If 
X -A , i >4 ? ' »V« U C 
X -0 # >/'iU7M? o r 

x - 0 . 

X — 0 , ViMvi 0 X 9 ? 

X -0. 

^OT A ( T,^L r 1 ( 1 > 
X ~n , -*.if,V*i) 0 4j r 
X 
X 
X 
x 
X 
X 
X 
X 
X 


— 0 • 4 i-l r «r> 1 1 x S- 
- 1 ' . 4 ?X *MI4;,r 


4 4 1 

i ‘>n‘, 

4 7 Mi *>P l 1 - 
t*J0:»-,P9q- 
iilAU ^ (1 t 

i9l 4-4 1 7 *' 

• VMI i?o./r 0 

( T 1 L r 1 ( I ' 

,nc,xc,.,r 9 

1 *F 

, ;7 .i*i7?.ir 
,ol ,>>XO;,r 
, ) 4 t | c » *? f > •* 

| 7 l X 4 4 / .; 

:j< > - 1 7 1; - 
,4 1 *1 •»*><! ) ■ 

,f,s7 'Om-- 
( T A * l r J ( l 1 
7 7 4 r % 4 • 1 / S’ M 

r > 7 **.>=; 0 


n 

o 

Ot 

10 

0 

•> 


-0. 

-a . 

-0 . 

-a. 

-n. 

-n . 

A t A 

x 

x - 0 . 

X -A. 

X -1. 

X 

x -O, 

X -M . 

x -n , 
x . 

X -f! . 

•XX. T A 
X -f* , 

X 

x . 7.10991 o 

X -0 , 7 | 1 j*i US.,.* O 
"AT A ( TA :i L r v I I . 

X -»• , 1 'n0». A 7 1 ‘ - 
X - 9 . 1 1 7 !. 0 «'S / ' 

x -<» , .» .144 s 7 n »*. 

X -0 , •«ti 1 f.7J 
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